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ABSTRACT 

We report the detection of a strong, organized magnetic field in the secondary component of 
the massive 08III/I+07.5V/III double-lined spectroscopic binary system HD 47129 (Plas- 
kett's star), in the context of the Magnetism in Massive Stars (MiMeS) survey. Eight inde- 
Ov pendent Stokes V observations were acquired using the ESPaDOnS spectropolarimeter at the 

Canada-France-Hawaii Telescope and the Narval spectropolarimeter at the Telescope Bernard 
Lyot. Using Least-Squares Deconvolution we obtain definite detections of signal in Stokes V 
in 3 observations. No significant signal is detected in the diagnostic null (AO spectra. The 
Zeeman signatures are broad and track the radial velocity of the secondary component; we 
therefore conclude that the rapidly-rotating secondary component is the magnetized star. Cor- 
rS . recting the polarized spectra for the line and continuum of the (sharp-lined) primary, we mea- 

■ sured the longitudinal magnetic field from each observation. The longitudinal field of the 

secondary is variable and exhibits extreme values of -810 + 150 G and +680 + 190 G, im- 
plying a minimum surface dipole polar strength of 2850 + 500 G. In contrast, we derive an 
upper limit (3<x) to the primary's surface magnetic field of 230 G. The combination of a strong 
magnetic field and rapid rotation leads us to conclude that the secondary hosts a centrifugal 
magnetosphere fed through a magnetically confined wind. We revisit the properties of the op- 
tical line profiles and X-ray emission - previously interpreted as a consequence of colliding 
stellar winds - in this context. We conclude that HD 47129 represents a heretofore unique 
stellar system - a close, massive binary with a rapidly rotating, magnetized component - that 
will be a rich target for further study. 
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Table 1. Log of observations and magnetic diagnosis of Plaskett's star. Phases correspond to the orbital ephemeris of lLinder et al. 1 2008), for which phase 0.0 
denotes the time of primary conjunction. In the "Det?" column, ND=No Detection, DD=Definite Detection. Detection probabilities in Stokes V (Py and null N 
Pn) are given in percent. Longitudinal magnetic fields (Be for Stokes V, N{ for the diagnostic null) correspond to t he LSD profiles of the secondary component, 
using the mean Stokes / profile of the secondary derived from the disentangled spectrum of lLinder et al. 12008) renormalised to the binary continuum. The 
longitudinal field detection significance (z = B;/<x) is given for both Stokes V (zb) and null N (zn). 
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1 INTRODUCTION 

HD 47129 (Plaskett's star) is an o ptically bright (V = 6.06), high- 
mass (M tot sin 3 i = 92.7 + 2.7 M^ilLinder et al.ll2008l) . short-period 
(P orb = 14.396257 ± 0.0 0095 d; llinder et al.ll2008l) . approximately 
circular dSticklandl 1987b non-eclipsing O+O SB2 system. It is con- 
sidered to be a probable member of Mon OB2, a large association 
located at a heliocentric distance of 1 .4 - 1 .7 kpc (Ch en et al J 20071) 
contai ning two or three subgroups of OB stars a nd stellar aggre- 
gates [Sing h & Naranarj|l979l : lli & Smith! |2005|) . These compo- 
nents are characterized by a relatively large range of character- 
istic ages, from 0.2-0.6 Myr (the dynamic al age of NGC 2244, 
the youngest component of the association; M athews) 1 1967t) to 2 
Myr (the main sequence turnoff age of NGC 2244; Park & Sung 
2002h to 4 Myr ( the dynamical age of the H i shell of the Rosette; 
Chen et alj2 007) to 20-25 Myr humeri 197(3) for the most evolved 
stars in the association. 

According to the analysis of llinder et al.l d2008l) . the sys- 
tem is composed of an 08III/I primary (HD 47129A1) of mass 
(45.4 ± 2.4) sin 3 i M and v sin ( = 60 - 75 km s 1 , and an 
07.5V/III secondary (HD 47129A2) of mass (47.3 ± 0.3) sin 3 ( M Q 
and v sin ( = 230 - 310 km s . Despite their very similar masses, 
the pri mary appears to be ~ 2x brighter than the secondary in the 
optical (Lin der et al.l 2008). The system inclination i = 71+9° 
was estimat ed by Rudy & Herman ( 1978) using linear polarimetry. 
iRudv & Herman! dl978l) and lBagnuolo et all dl992b noted that an 
inclination significantly larger than ~ 70° would result in eclipses 
(which are not observed), while an inclination significantly smaller 
would imply improbably large masses of the components. There- 
fore, with a total mass of approximately 1 12 M Q , HD 47129 is one 
of the most massive known O-type binaries. 

The system is clearly chemically peculiar: according to 
Linder lltall d2008l) . the primary is strongly N enhanced (16 times 
solar) and C depleted (3% of solar), while the secondary is N de- 
pleted (20% of solar) and He overabundant (1.5 times solar). The 
peculiar chemistry, in combination with the mass/luminosity mis- 
match and rapid rotation of the secondary, has led investigators 



dBagnuolo et al.lll992l ; iLinder et alj|2008b to speculate that Plas- 
kett's star is a post Roche lobe overflow (RLOF) system. 

The system was first identified as an X-ra y emitter by 
the RO SAT All-Sky Survey dBerghoefer et al.ll 19961) . [Linder et all 
(2006) reported analysis of XMM-Newton observations of Plas- 
kett's star, revealing the system to be a hard, luminous and variable 
X-ray emitter (kT max ^ 1.4 keV, logL x /L bol = - 6.0). These ob- 
serva tions support the historical interpretation (e.g. Iwiggs & Giesl 
1 1992b of the system as a colliding-wind binary (CWB). 

HD 47129 was observed in the context of the Magnetism in 
Massive Stars (MiMeS) survey of magnetic fields in Galactic OB 
stars. This target was flagged as a probable Zeeman detection fol- 
lowing preliminary analysis via the MiMeS pipeline. The primary 
goal of this paper is to report the firm detection of a strong magnetic 
field in the secondary component of the system. This discovery 
leads directly to a re-examination of the optical and X-ray prop- 
erties of Plaskett's star as a Magnetically-Confined Wind Shock 
(MCWS) system or CWB/MCWS hybrid. The latter classification 
would make HD 47129 unique among the known massive magnetic 
stars. 



2 OBSERVATIONS 

Fifteen high resolution (R ^ 65 000) spectropolarimetric (Stokes / 
and V) observations of HD 47129 were collected with ESPaDOnS 
at the Canada-France-Hawaii Telescope (CFHT) and Narval at 
the Telescope Bernard Lyot (TBL) between 3 Feb and 25 Mar 
2012. Each spectropolarimetric sequence consisted of four indi- 
vidual subexposures, each of 600 s duration (for ESPaDOnS) or 
1 200 s duration (for Narval), taken in different configurations of 
the polarimeter retarders. From each set of four subexposures we 
derived Stokes / and Stokes V spectra in the wavelength range 

3 670 - 10 000 A f ollowing the double-ratio procedure described by 
iDonati et all dl997l) . ensuring in particular that all spurious signa- 
tures are removed to first order. Diagnostic null polarization spectra 
(labeled AO were calculated by combining the four subexposures in 
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Figure 1. Stokes / (bottom), V (top) and diagnostic null N (middle) LSD profiles corresponding to the highest SNR polarized spectra (From left: JDs 2455961 
(DD), 5966 (ND), 5967 (DD), 5969 (DD) and 6012 (ND)). No significant signal is detected in N. V and N have been scaled (by 300x) and shifted for display 
purposes. Note the complexity and variability of Stokes /. 



such a way that polarization cancels out, allowing us to verify that 
no spurious signals are present in the data (for more details on the 
definition of N, see lDonati et alJll997T) . All frames we re processed 
using the automated reduction package Libre ESpRIT (Don ati et al.l 
11997?) . 

All spectra acquired during a single night of observation were 
co-added following reduction to increase the signal-to-noise ratio 
(SNR). This ultimately resulted in 8 independent spectropolarimet- 
ric observations (i.e. 4 co-added ESPaDOnS spectra and 4 individ- 
ual Narval spectra). The peak SNRs per 1.8 km s pixel in these 
spectra range from 600-3000. The log of observations is reported 
in Tabled 



3 DETECTION OF A MAGNETIC FIELD IN 
PLASKETT'S STAR 

As described by e.g. iLinder et~al", I J2008I) . the optical spectrum of 
HD 47129 is extremely complex. Due to the broad lines of the 
secondary, spectral lines of the two components are never well- 
separated in velocity during the +200 km s amplitude radial ve- 
locity (RV) orbits of the two stars. Many lines exhibit multiple, 
variable (in position and profile shape) emission and absorption 
components. These spectral characteristics make a detailed diag- 
nosis and interpretation of the magnetic properties of the system 
very challenging. In this first analysis, we employ Least-Squares 
Deconvolution (LSD; Don ati et al .1 19971) and take adv antage of the 
disent angled spectra of the components computed bv lLinder et al.l 
(2008) to aid in the interpretation of our results. 

Examination of the best-quality co-added Stokes V spectra re- 
veals no strong evidence of Zeeman signatures in individual spec- 
tral lines of the primary or secondary (although, as described above, 
the complexity of the spectrum impedes such an identification). 
To improve our detection capabilities, LSD was applied to the re- 
duced spectra. To begin, we adopted a line mask based on the LTE 
synthetic spectrum of an 08 giant. As in the case of other mag- 
netic O-type stars with complex spectra (e.g. iDonati et alj|2006l : 
iMartins et al.l201oh . we found that a relatively small subset of lines 
(12) contribute usefully to the detection of a Zeeman signature. We 
ultimately arrived at an optimal line mask, yielding the most signifi- 
ca nt Zeeman detection s, with content similar to the mask employed 
bv lDonatietailfeOOfj) for HD 191612. For more detail regarding 
line mask optimization for O-t ype stars wi t h complex spectra, the 
reader is invited to consult, e.g. Wade et alj \20 111 12013) . 

Extraction of LSD profiles from the 8 spectropolarimetric ob- 



servations using the adopted mask yields 3 definite detections and 
5 non-detections of significant signal in the St okes V LSD profile s 
(according to the detection criteria described bv lDonati et al.ll997h . 
In no case is any significant signal detected in the diagnostic null. 
The LSD profiles corresponding to the five highest-SNR spectra - 
i.e. 3 definite detections as well as two higher-SNR non-detections 
- extracted on a velocity grid spanning 2000 km s _I sampled at 
27 km s _1 , are illustrated in Fig.Q] Note that while a majority of 
our observations correspond formally to non-detections, the Stokes 
V profiles of most observations are visibly more structured (and ex- 
hibit larger detection probabilities) than the corresponding A' pro- 
files (see, for example, the second and fifth panels of Fig. [T] as well 
as Table 1). 

A first examination of the LSD profiles in Fig.Q~]reveals that 
the Stokes V Zeeman signatures are generally much broader than 
the sharp, easily identified Stokes / profile of the primary star. This 
strongly suggests that the Stokes V profile has its origin in the spec- 
trum of the (broad-lined) secondary star. Unfortunately, while the 
primary's contribution to the LSD Stokes / profiles is distinctive, 
the position and range of the broad secondary profile is more diffi- 
cult to determine. To aid us in this task, we used the disentangled 
spectr a of the primary and secondary computed by ILinder et al.l 
(2008). For each of the 8 observed orbital phases, we applied ap- 
propriate RV shifts to the disentangled s pectra corresponding to 
the orbital solution of ILinder et alj J2008I) . We then extracted the 
Stokes / LSD profiles from these shifted spectra in the same man- 
ner as was applied to the ESPaDOnS and Narval Stokes / spectra. 

In principle, these LSD profiles should be representative of 
the individual components' LSD profiles at the observed orbital 
phases, assuming minimal intrinsic variability of the stars. To check 
the limitations of this assumption, we computed the binary LSD 
profiles expected at the observed phases by weighting the individ- 
ual profiles according to the optical brightness ratio (L pY JL xc = 
1.9; lLinderetal 1 12008), summing them and renormalising them to 
the combined binary continuum. Comparing these to the observed 
binary LSD profiles (from the ESPaDOnS/Narval spectra; Fig. [2j 
we find an acceptable agreement. We therefore conclude that while 
the primary and secondary are variable, the magnitude of this vari- 
ability is not sufficient to introduce fundamental uncertainties into 
the following analysi s. This conclusion is consistent with that of 
iBagnuolo et al. I dl992l) . who report that "the secondary contribution 
to the cross-correlation profiles is relatively constant with orbital 
phase." 

Comparing the observed Stokes V profiles with the position 
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Figure 2. Observed and modelled Stokes / LSD profiles of HD 47129 (in chronological order, left to right and top to bottom). Top row in each frame: LSD 
profile of the secondary shifted to the predicted secondary RV at the time of observation. Middle row in each frame: LSD profile of the primary shifted to the 
predicted primary RV at the time of observation. Bottom row in each frame, thick (Mack) line: Observed LSD profile of HD 47129. Bottom row, thin (red) 
line: model LSD profile obtained from weighted and renormalised addition of the LSD profiles of the individual components extracted from the disentangled 
spectra. Orbital phase is indicated in the lower left corner of each frame. 



and width of the primary's mean Stokes / profile confirms our 
initial assessment: the primary spectrum is not the origin of the 
detected Zeeman signatures. On the other hand, the position and 
breadth of the secondary's mean Stokes / profile agree well with 
the Stokes V profiles at all phases. The data therefore support the 
hypothesis that the secondary component of HD 47129 is a rapidly 
rotating, magnetized O-type star. 



4 MAGNETIC FIELD CHARACTERISTICS 
4.1 Measurement of the longitudinal field 

Having concluded that the secondary component is responsible for 
the observed Stokes V Zeeman signature, we wish to constrain the 
magnetic strength and geometry of the star at the level of the photo- 
sphere (where the lines employed in the LSD mask are principally 
formed). In the following we assume that the primary star con- 
tributes negligibly to the observed circular polarization. While this 
assumption is consistent with the observations, we cannot strongly 
test it given that the primary's profile is always blended with that of 
the secondary. On the other hand, it does not appear that the primary 
contributes to Stokes V at a level comparable to the secondary. 

The simplest quantitative measure of the stellar magnetic field 
from circular polarization measurements is the mean longitudinal 
magnetic field (or simply longitudinal field), computed from the 
first-order m oment of the Stokes V profile about it s velocity centre- 
of-gravity v jDonati et alJl997l ; IWade et alj2000l) according to the 
equation: 



2.14 x 10 1 

Xzc 



(v - v )V(v) dv 



J[l-/(v 



(1) 



)]dv 



In Eq. (1) V(v) and /(v) are respectively the continuum- 
normalised Stokes / and V profiles as a function of velocity, v. 
The wavelength A is expressed in nm and the longitudinal field Be 
is in gauss. The wavelength and Lande factor z correspond to the 
weighting factors used in calculation of the LSD profiles (500 nm 
and 1.2, respectively). In our mask, the mean SNR- weighted wave- 
length and Lande factor were 511 nm and 1.11, respectively. 

Note that the integral in the denominator of Eq. (1) represents 
the equivalent width of the Stokes / profile (in velocity units). To 
infer the longitudinal magnetic field of the secondary as a single 
star, we therefore require the secondary's line profile independent 
of the primary's. We attempted to subtract the primary profile from 
the combined profile at each phase using the primary's LSD profile 
from the disentangled spectrum. However, the primary's line is suf- 
ficiently sharp and deep that relatively small variability or errors in 
RV seriously distort the resultant secondary profile. Ultimately, we 
decided to replace the binary Stokes / profile with the secondary's 
LSD profile from the disentangled spectrum. The accuracy of this 
choice is limited in two ways. First, the secondary profile is vari- 
able l lLinder et alJl2008T) : it is therefore likely that its equivalent 
width also varies. Secondly, the profile shows important emission 
in its wings. In all likelihood this emission is formed in the (proba- 
bly magnetically-confined) wind, above the photosphere of the star, 
where the magnetic field is much weaker. Consequently, we assume 
that the emission does not contribute significantly to the Stokes V 
profile (this assumption is revisited in Sect. 7). On the other hand, 
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Figure 3. Stokes / (bottom), V (top) and diagnostic null N (middle) LSD profiles corresponding to selected observations (as in Fig.[TJ used for the calculation 
of the longitudinal field of HD 47129. V and N have been scaled and shifted for display purposes. The ob served bin ary Stokes / profile has been replaced by a 
shifted and renormalised LSD profile of the secondary star, computed using the disentangled spectrum of lLinder et alj 420081) . The integration bounds used to 
evaluate Eq. (1) are illustrated with vertical dotted lines. 



it reduces the equivalent width of Stokes / if included in the range 
of velocity used to evaluate Eq. (1). 

Therefore to evaluate Eq. (1) we used the LSD profile from 
the disentangled secondary spectrum, shifted to the RV appropriate 
to each observation, and renormalised to th e continuum of the real 
system (using the optical brightness ratio o f lLinder et aljioolh (Al- 
ternatively, we could have renormalised the Stokes V profile to the 
continuum of the secondary. It can be shown that for the calculation 
of Be these operations are equivalent.). Eq. (1) was evaluated in the 
velocity range ±300 km s _1 symmetric about vo. The velocity range 
was selected based on the reported v sin i of the secondary, as well 
as visual examination of the Stokes / and V profiles. We find that 
this range corresponds well to the central absorption region of the 
line profile. We recognize that emission may partially infill the ab- 
sorption line, but we have no reasonable way to estimate the impor- 
tance of such infilling. However, if our hypothesis that the emission 
results from a rigidly-rotating centrifugal magnetosphere is correct, 
then emission will be confined to velocities larger than v = ±v sin i, 
and will therefore not significantly infill the line profile. A selection 
of the LSD profile sets used for determination of the longitudinal 
field are illustrated in Fig. [3] In this figure we show the Stokes V 
and null N profiles (identical to those in Fig.QJ, and the shifted and 
renormalised Stokes / LSD profile of the secondary derived from 
the disentangled spectrum. The integration bounds used to evaluate 
Eq. (1) are illustrated with vertical dotted lines. 

A useful product of the calculation of the longitudinal mag- 
netic field is the determination of the centre-of-gravity vo of Stokes 

V at each phase. These measurements provide a quantitative veri- 
fication of our tentative conclusion that the magnetic field is asso- 
ciated with the secondary component. As illustrated in Fig. [4] the 
vo measurements track the RV of the secondary faithfully. In fact, 
these measurements provide an independent verification of the RV 
variation of the secondary that is independent of the disentangling 
procedures employed in past studies. Clearly they do not track the 
primary's RV curve well at all. 

It is therefore clear that the dominant contribution to Stokes 

V comes from the secondary. Nevertheless, a weak contribution 
to Stokes V - even above the level of the noise - could possibly 
come from the primary, and remain hidden in the stronger Stokes 

V profile of the secondary. Due to the primary's sharper spectral 
lines, a magnetic field is intrinsically more easily detectable in this 
star. Based on the ratio of projected rotational velocities of the two 
stars, we would roughly expect the primary's Stokes V signature to 



have an amplitude equal to that of the secondary if its longitudi- 
nal magnetic field were one-fourth as large as the secondary's (i.e. 
(vsin;)2/(Vsin;)i = 70/305 - 0.25). Based on our analysis, we 
consider such a signal to be ruled out at 3cr confidence. Taking into 
account the optical brightness ratio of the two stars, this yields a 3cr 
upper limit on the longitudinal field of the primary on JD 2455969 
(i.e. when the secondary's field was detected with the greatest sig- 
nificance) of about 70 G. 



4.2 Variability of the longitudinal field 

The longitudinal field measurements of the secondary are reported 
in Table 1. They are observed to vary from an extreme positive 
value of +680± 190 G to an extreme negative value of - 810± 150 G. 
The variation is significant at about 6.2<x, and the change of the sign 
is reflected in the changing shape of the Stokes V profile. The lon- 
gitudinal field measured on consecutive or near-consecutive nights 
- in particular on 9/10/12 Feb - is observed to be significantly dif- 
ferent, indicating that the variability timescale is rather short, likely 
of order 1 d. 

We have examined the period content of the longitudinal field 
measurement s using the FFT and clean procedures as implemented 
in CLEAN- NG (Gutie rrez-Soto et alj |2009). The only significant pe- 
riod present in the 8 longitudinal field measurements corresponds 
to 3.6 ±0.1 d, and yields an approximately sinusoidal variation 
of B ( , with a reduced x 2 of 1-37. (In contrast, folding the data 
with the orbital period yields a substantially worse fit, with reduced 
X 1 = 7.46.) However, given the relatively low SNR of the variation 
and the small number of data points, the reliability of this period is 
unclear. However, it is still instructive to compare it to the expected 
rotational periods of the stars. 

Rotational periods of the components have been computed us- 
ing raduinferred from the temperatures and luminosities reported 
by iLinder et alj J2008I). The uncertainties on the luminosities de- 
rived bv lLinder et al.l d2008l ) do not account for the uncertain dis- 
tance to Mon OB2. To account for the distance modulus uncer- 
tainty, we adopt a luminosity uncertainty 0.1 dex. We note that this 
calculation does not take into account any rotational effects (e.g. 
oblateness or gravity darkening), nor does it include any mutual 
gravitational defor mation of the stars. The projected rotational ve- 
locities derived bv lLinder etail d2008l) using the Fourier technique 
on He i and He n lines range from 60 - 75 km s for the primary 
and 230-310 km s _1 for the secondary. For the purposes of our cal- 
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culations, for t he primary we adop t vsirn = 70 ± 10 km s -1 based 
on the results of Linder et al. ( 20Qa), but w eighted to slightl y highe r 
velocity taking into account the result of iBagnuolo et al ] l ll992h . 
For the secondary, we adopt 305 ±15 km s _1 . Here we have take n 
into account only the results for He n lines bv llinder et alj i2008t) . 
since the He i lines are visibly affected by emission and result in 
substantially lower vsini values. Th e adopted vsin ; of th e sec- 
ondary is consistent with the result of Bagn uolo et alj dl992h . 

For the convenience of the reader, the relevant data are sum- 
marized in Table 2. We derive the radii of the components: 14.1 ± 
3.3 R e for the primary, 10.7 ± 2.6 R a for the secondary. We obtain 
rotational periods P rot ,i / sin i = 10.2 ± 4.0 d and P tas j,l sin i = 1 .8 ± 
0.5 d. If we assume that the rotational and orbital angular momenta 
are a lignecfl and adopt an inclination i = 71° ( Ru dy & Hermanl 
Il978h , then P mtA = 9.6 ± 3.7 d and P nu2 = 1.7 ± 0.5 d. The 
longest possible secondary rotational period consistent with the de- 
rived radii and published v sin i (i.e. for i = 90°)is 1.8 + 0.5 = 2.3 d. 
Therefore the 3.6 d period inferred from the FFT is not straightfor- 
war dly compatib l e with rotational modulation by the secondary. 

iMahv et alj < T201ll) analyzed 34 days of continuous CoRoT 
photometry of Plaskett's star. The two most significant periods de- 
tected correspond to the orbital period of 14.4 d, and 1.22 d. They 
report no immediate explanation for the 1.22 d period (although we 
do note that it is approximately one-third of 3.6 d), but conclude 
that this is not likely to be related to the rotational period of either 
star. This could only be the base if the stellar rotation axis were 
inclined to our line-of-sight by 43°, in which case this value would 
be consistent with the inferred radius and v sin ( of the secondary. 
On the other hand, as discussed in Sect. 1, if the rotation and or- 
bital axes are aligned, such an inclination would imply masses of 
th e components signif icantly (2.6 times) larger than those adopted 
bv lLinder etal] fe008l) . 

Two a dditional periods, 2.51 d and 2.71 d, are reported by 
iMahvetalJbOllh . These are similar t o the period of varia tion of 
the Her wings of P = 2.86 d reported bv lWiggs&Giel ( ll992h . and 
somewhat longer than the longest admissible secondary rotation pe- 
riod. Otherwise we are not yet able to offer any useful interpretation 
of them. 

Given the relatively large amplitude of the longitudinal mag- 
netic field of HD 47129A2, we are confident that further spectropo- 
larimetric observations will yield a direct determination of the sec- 
ondary's rotation period. 

4.3 Magnetic strength and geometry 

Although the available magnetic data are insufficient to develop 
a unique model of the magnetic geometry of HD 47129A2, it is 
reasonable to assume (similar to other magnetic O stars) that HD 
47129A2 is an oblique rotator, and (given the large value of the 
longitudinal field and its strong variability) that the magnetic field 
has a relatively simple structure - i.e. it is domin ated by t he dip ole 
component. In this case we can use Eq. (1) of IPreston d 19671) to 
estimate the minimum possible polar strength of the surface dipole. 
Using a limb-darkening coefficient of 0.3 and the peak longitudinal 
field l-Bflmax = 810 ± 150 G, we compute a minimum polar field 
strength of the dipole B d = 2850 ± 500 G. 

While we are not able to unambiguously determine the period 

1 While this appeal's to be a reasonable assumption for HD 47 1 29 (see Sect. 
7), other close, high-mass bina ries exist for which this may not be the case 
IVillar-Sbaffi et alj2005ll200^ . 



Table 2. Summary of physical properties of HD 47129 used to compute 
component rotation al periods. T e g, log L„ and v sin i are values reported by 
iLinder et alj |200§t) (except that the luminosity uncertainties have been in- 
creased from 0.04 to 0.1 dex, as described in the text). R„ and P rot are cal- 
culated from those values. 
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Orbital Phase 

Figure 4. Triangles: Radial velocities of the primary component measured 
from ESPaDOnS and Narval Stokes / LSD profiles. Circles: Centres-of- 
gravity of the Stokes V LSD profiles. The dashed and solid lines repre- 
sent respectively the theoretical RV variations of the pri mary and s econdary 
components according to the orbital model of lLinder et alj (2008). Clearly 
the Stokes V profiles follow the secondary's RV variation, further establish- 
ing the association of the magnetic field with the secondary component. 

of variation of the longitudinal field, it is clearly detected to vary 
and to change sign. This indicates that both magnetic hemispheres 
are visible during the star's rotation. This implies that at least one 
of the two angles specifying the magnetic field geometry (the incli- 
nation i of the rotation axis to our line-of-sight, and the obliquity /? 
of the magnetic axis relative to the rotation axis) is relatively large, 
and that i+p> 90°. 

Applying similar reasoning to our upper limit on the longi- 
tudinal field of the primary, we obtained a 3cr upper limit on the 
primary's surface dipole component of approximately 250 G, i.e. 
more than lOx smaller than that detected in the secondary. 



5 THE WIND-MAGNETIC FIELD AND WIND/WIND 
INTERACTIONS 

As discussed by lud-Doula & Ow ocki (2002), the global compe- 
tition between a magnetic field and stellar wind can be charac- 
terized by the so-called wind magnetic confinement parameter 
77, = P>l q Rl/ Mv m , which depends on the star's equatorial surface 
field strength (S cq = 5d/2), radius (/?,), and wind terminal momen- 
tum (Mv„). For a dipolar field, one can identify an Alfven radius 
^Aif — ll 14 R*, representing the radial extent of strong magnetic 
confinement of the wind. In the equatorial plane, above /?Aif, the 
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wind dominates and stretches open all field lines. But below R Ali , 
the wind material is trapped by closed field line loops. In the case 
of a rapidly rotating star like HD 47129A2, centrifugal support re- 
sults in a long residence time of confined plasma located above the 
Kepler co-rotation radius R^ tp = (GM/w 2 ) 1 ' 3 , where u> is the rota- 
tional frequency of the star (i.e. a> = v Iot /R,). 

Although the primary's mass los s rate and wind term inal speed 
were determined observationally by Linderet al. (2008), the sec- 
ondary's were not. In fact. lLinder et alj d2008l) remark on the poor 
match of their CMFGEN spectra to optical and UV spectra, and the 
puzzling absence of the classical wind lines in the secondary's UV 
spectrum. 

To estimate the wind momentum of the secondary, we used 
the adopted parameters given in Table [2] and masses of 54 M Q 
for the primary and 56 M Q for the secondary (i.e. the projected 
masses cited in Sect. 1 evaluated for an orbital inclination of 71°). 
While outside the scope of this paper, we note that the spectro- 
scopic masses determined using the inferred temperature and lumi- 
nosity are significantly (a factor of two for the primary, and 4 for 
the secondary) smaller than those determined dynamically. This is 
a serious discrepancy that may be tied both the the peculiar spectra 
and to the exotic history of the system, discussed further in Sect. 7. 
Whatever its origin, this discrepancy demands urgent attention. 

We determined the the oretical wind term inal velocity and the 
mass-loss rate according to I Vink et"al . (2000,l200l|). We obtained 
vv2 = 3500 km s" 1 and M 2 = 7.1 x lO^Moyr 1 . Although the 
calculated terminal speed is higher than that m easured for the pri- 
mary (Vco,] = 2450 km s -1 ; ILinder eTa l. 2008), t he m ass loss rate 

20081) . (For com- 
2980 km 
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is lower (Mi = 1.7 x 10~ 7 M G yr~ 
parison, similar calculations for the primary yield v, 
s -1 and Mi = 2.9 x 10~ 7 M Q yr~', which are in reasonable agree- 
ment with the measured values.) From a physical perspective, M of 
the secondary is expected to be lower because the star is less lumi- 
nous than the primary, but v m is larger because the escape velocity 
(determined by the surface gravity, GM,/Rl) is larger. 

These wind characteristics combined with the estimate of 
the surface magnetic field strength imply n, = 680: the mag- 
netic field strongly confines the wind, dominating it dynamically 
to R A \f = 5.4 in the magnetic equatorial plane. If we underes- 
timated the wind momentum of the secondary by a factor of e.g. 
two, this would result in an overestimation of R Alf by 1.3 R,. We 
conclude that the uncertainty in R A \f is dominated by the system- 
atics in the wind momentum calculation, and not the formal errors 
of the stellar parameters. In addition, we recall that Sj represents a 
lower limit to the magnetic field strength, and is influenced by any 
systematic uncertainty in the optical brightness ratio. 

HD 47129A2 is the most rapidly-rotating known magnetic O- 
type star. We therefore expect rotation to have an important role in 
determining the structure of its magnetosphere. Using the values 
in Tableland adopting sin; = 1 (i.e. assuming the rotation is the 
slowest permitted by v sin f) we compute the Kepler co-rotation ra- 
dius of the secondary, Rk cp < 2.2 R,. That is already below our best 
estimate of the lower limit on R Alf . This implies that HD 47129A2 
hosts a region, between R ^ 2.2 R, and R = 5.4 R, (i.e. from 
1.2-4.4 R, above the stellar surface) in which the wind is both mag- 
netically confined and centrifugally supported. Thus we conclude 
that the secondary hosts a centrifugal magnetosphere - the first such 
configuration amongst the known magnetic O-type stars. 

Well-known examples of centrifugal magn etospheres exist 
amongst the magnetic B -type stars (e.g. or Ori E; iTownsend et al.1 
120051 : loksala et ai]l2012h . The prototypical optical spectral signa- 
ture of these configurations is emission in B aimer lines at high 



velocities (i.e., several times vsini, as a consequence of magneto- 
spheric plasma maintained in rigid rotation with the star, high above 
the stellar surface). In the case of <x Ori E, for examp le, the emis- 
sion e xtends to ±600 km s _I , or about ±3.75vsin; dOksala et al.l 
l2012h . (For reference, the Alfven and Kepler radii for cr Ori E are 
estimated to be 1.9 R, and 31.1 R,; Petit et al., in preparation.) In 
the case of HD 47129A2, analogous emission at high velocities is 
observed in many spectral lines (pr incipally lines of H and He 1; see 
Figs. 1 and 8 of iLinder et alj[20o3) . Fig[5]shows the LSD profile 
extracted from the disentangled spectrum of the secondary compo- 
nent, as well as the He 1 /14471 line from the disentangled spectrum. 
Vertical dotted lines indicate ±v sin i (= ±305 km s -1 ) and the es- 
timated limits of the emission (at -900 km s _1 and +710 km s ). 
The emission therefore extends to approximately 2.3 to 2.9 times 
v sin (. These values are certainly compatible with the Alfven and 
Kepler radii calculated above, and with the assumption (from Sect. 
4.1) that the magnetic field in the emitting material contributes neg- 
ligibly to the observed Stokes V profil^]. 

Attribution of this emission to a magnetosphere, rather than 
stellar wind collision, makes sense given th at no simila r emis- 
sion is associated with the primary. In fact, ILinder et alj fe008h 
propose based on Doppler tomography that the secondary is sur- 
rounded by a rotationally-flattened wind, concentrated in the equa- 
torial plane. They comment that this could explain the apparently 
stronger wind of the secondary (based on emission line strengths 
and X-ray absorption). It seems reasonable, given the detection of 
a magnetic field, to identify an equatorial enhancement of the sec- 
ondary's wind with plasma confined by the star's magnetosphere. If 
the secondary's wind is structured into a magnetosphere, this may 
also help to explain the highly peculiar UV wind lines, in an alogy 
with other magnetic early-type stars (e.g. lHenrichs et aU2005l) . 

The wind momentum ratio M 2 v , oo,2/ / Mi l Vi of the components 
is calculated to be 0.6 (using the measured wind properties of the 
primary) or 0.3 (using the computed wind properties of the pri- 
mary). In both cases, we conclude that the primary's wind domi- 
nates. This implies that the primary's wind wraps around the sec- 
ondary component, and places the wind collision at a distance of 
3.1 R, (or 1.6 R,) from the secondary. Thus the classical wind col- 
lision, in the absence of the dynamical influence of the magnetic 
field, is predicted to occur inside the Alfven radius of the secondary 
component. It could well be the primary's wind is actually deflected 
by the magnetic field of the secondary (just as the same magnetic 
field limits the wind of the secondary), in analogy to the deflection 
of the solar wind around planetary magnetospheres in the solar sys- 
tem. Or the primary's wind might be channelled along open field 
lines, effectively feeding the secondary's magnetosphere. In both 
scenarios the winds never reach a classical CWB collision, at least 
in a direction close to the line joining the two stars. Based on this 
reasoning, we should ask ourselves if the CWB scenario is in fact 
a good description of the HD 47129 system. 

It is also interesting to consider that the magnetic field of 
the secondary is likely to be tilted relative to the orbit axis, and 
therefore that the magnetospheric structure will periodically rotate 
around or through the orbital plane, potentially generating a highly 
dynamic interaction with the wind of the primary. When the sec- 
ondary's magnetic pole is directed at the primary, the relatively 
unimpeded wind of the secondary could produce a more typical 



2 Assuming a dipole field falling off as 1/r 3 (as will be the case inside the 
Alfven radius), the field at these distances is 12-24 times weaker than at the 
stellar surface. 
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Table 3. Fit results for vm Lya. Uncertainties are reported for 95.3% 
confidence interval. 
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To the extent that it can be constrained, the broadening is compara- 
ble to that observed for O vm Lyo- in the spectrum of 9 [ Ori C. 

Thus, we find that both CWB and MCWS models for X-ray 
emission are consistent with the observed X-ray properties of Plas- 
kett's star. 



Figure 5. He i A441 1 line (thick black line) and LSD profile (thin red line) 
extracted from the disentangled spectrum of the secondary star. Vertical 
lines indicate ±v sin i and the apparent extent of the wing emission. 



colliding wind scenario. But when the secondary presents its mag- 
netic equatorial region to the primary, this should result in closed 
magnetic loops (and confined wind plasma of the secondary) being 
impacted by the primary wind. These differences could potentially 
influence the X-ray light curve. Basically, one might expect mod- 
ulation according to the rotation period of the secondary, with two 
peaks and troughs per rotation cycle, depending on geometry of the 
magnetic field (i.e. the angles i and /3). 



6 X-RAYS 

The properties of t he X-ray emission of Plaskett's star reported by 
iLinder et al.l J2006I) based on their observations with XMM-Newton 
have been interpreted as consistent with a CWB scenario, which 
was reasonable given the knowledge of the system at that time. 

However, the relatively high X-ray luminosity (log L x /L bo ] 6.0) 

an d temperature (T ma x ~ 2.5 keV) are comparable to those of 9' Ori 
C dGagne et alj2005l) and other known magnetic O stars, and could 
be consistent with either MCWS or CWB scenarios. 

We reana lyzed the archival XMM-Newton RGS spectra of 
Plaskett's star dlinder et alfeoOfj) to obtain constraints on two key 
parameters which could indicate the mechanism of its X-ray emis- 
sion. The ratio of the forbidden to intercombination lines (///) in 
He-like triplets can c onstrain the proxi mity of the X-ray emitting 
plasma to the stars teahnetalj|200lh . Doppler broadening and 
shifts of sp ectral lines can also provide contraints on the emitting 
geometry dOwocki & CoherfeOOll : [Henley et alj|2003l) . 

We reprocessed the XMM-Newton RGS data using SAS 
12.0.0. We used rgsproc to produce spectra and response matri- 
ces, filtering intervals of high background. 

We fit the Si xm, Mg xi, Ne ix, and O vn He-like triplets with 
a three Gaussian model. We find that the counting statistics of the 
data set are insufficient to provide significant constraints on the //( 
ratio of any of these species. 

We also fit several of the strongest lines with Gaussians to 
place constraints on the Doppler broadening and shifts. The best 
constrained line is O vm Lya, with best fit values given in Table[3] 



7 DISCUSSION AND CONCLUSION 

In this paper we have reported the firm detection of Zeeman signa- 
tures in high resolution Stokes V spectra of the massive binary HD 
47129 (Plaskett's star). We have demonstrated that the detected sig- 
natures are formed in the spectral lines of the rapidly rotating sec- 
ondary component. No comparable field is observed in the sharper- 
lined primary. Based on the large inferred value of the longitudinal 
magnetic field, and its significant variation, we conclude that the 
field is highly organized over the star's surface. The surface po- 
lar strength of the dipole component of the field must be at least 
2850 + 500 G. The field is observed to vary on timescales of days, 
compatible with rotational modulation of an oblique, static field 
in the rotat ing frame o f the secondary star (i.e. the oblique rotator 
model of Stibbs 1950). Although the data do not yet allow us to 
infer the period of variation of the longitudinal field, we find that 
the measurements do not phase coherently with the system's orbital 
period, nor with other periods derived from photometric or spectral 
variability. 

Computation of the magnetic wind confinement parameter 77, 
using theoretical values of the secondary's mass loss rate and wind 
terminal velocity indicate that magnetic channelling should confine 
the wind of the secondary to a distance of 5.4 R,. Coupled with 
the observed rapid rotation, this confinement is expected to pro- 
duce a long-lived centrifugal magnetosphere of stellar wind plasma 
in the region 2.2 - 5.4 R,. The high- velocity emission wings of 
the secondary's H and He 1 line profiles are quantitatively consis- 
tent with this scen ario, while the highly peculiar UV wind profiles 
(Linder et "al] |2008h of the secondary are qualitativel y consistent. 
The X -ray spectrum observed with XMM-Newton jLinder et al.l 
l2006t) is also consistent with the expected equatorial shocks due to 
collision of the wind flows channelled by the magnetic field from 
opposite magnetic hemispheres. This suggests that previous inter- 
pretations of Plaskett's star as a classical CWB may be too sim- 
plistic. Indeed, th e Doppler tom ographic maps of He 11 /14686 and 
Ho- computed by [Linder et alj J2008t> are quite unusual for collid- 
ing wind systems in the sense that they reveal a ring of material 
in velocity space. Such ring-like structures are often seen in Algol- 
type binaries (and there they are attributed to the accretion disk). 
Since apparently neither of the stars in HD 47129 (currently) fills 
its Roche lobe, the most likely explanation (in view of the detec- 
tion of the magnetic field) is that we are seeing the emission from 
the magnetically confined wind. The line profile variations are also 
quite atypical for CWBs: most colliding wind systems where the 
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bulk of the emission arises in the wind-wind interaction zone dis- 
play broad lines with changing profiles. In Plaskett's star, the broad 
component is relatively constant. What varies from one observa- 
tion to the next is a narrow component that moves with some phase 
shift compared with the primary star. It is this narrow component 
that dominates the Doppler map; this is because Fourier filtering (as 
applied to the spectra in the Fourier filtered back-projection method 
used to build the Doppler map), actually removes the strictly con- 
stant part of the spectrum, i.e. in this case part of the confined wind 
emission). All of this suggests that what we are seeing in this sys- 
tem is more complex than a classical wind-wind collision. 

As mentioned in Sect. 1, Plaskett's st ar is proposed to be a 
post mass transfer (i.e. post RLOF) system iBagnuol o et al.|[l~992l ; 
Linde r et al. 2008). The RLOF scenario is invoked to explain, at 
least qualitatively, the mass-luminosity mismatch of the system, the 
peculiar surface chemistry of the components, and the secondary's 
rapid rotation. In this model, the primary was initially the more 
massive star. As it evolved, it filled its Roche lobe. The secondary 
gained rotational angular momentum during the mass transfer pro- 
cess, likely via an accretion disk located in the plane of the orbit. 
The secondary would only have needed to accrete a f ew percent of 
its own mass to reach critical rotation^] (Packet 1981). During this 
event the angular momentum of the secondary would have been 
strongly modified, and likely forced into alignment with the axis of 
the orbit. Rapid circularization p robably also occurred during the 
RLOF event (e.g. lHurlev et alj|2002h . 

The RLOF mass transfer process would certainly have rep- 
resented a violent and transformative event in the life of the sec- 
ondary component. If we imagine that the magnetic field we ob- 
serve is a fossil field that existed prior to the RLOF event, we ex- 
pect that mass transfer would have resulted in strong shear in the 
outer envelope, seriously distorting the structure of the field in the 
outer envelope. The shear would have acted on the poloidal com- 
ponent, generating toroidal field. If the toroidal field became suf- 
ficiently strong, and if the shear was not frozen by back-reaction 
of the field, a Tayler instability would likely have resulted. At 
the same time, Maxwell stresses and the Lorentz torque would 
have acted to damp differential rotatio n along field lines (see 
iMestel et alJl 19881 ; iMathis & Zahnll2005l) while the instability sat- 
urated. Moreover, if the instability was was not able to regenerate 
both components of the initial field, the field would ha ve decayed 
by oh mic diffusion. Such processes are observed by IZahn et all 
J2007I). Therefore, if a fossil field was present at the be ginning (e.g. 
iBraithwaite & Nordlun J I200I ; iDuez & MathJsl l2010h . its charac- 
teristics would have been strongly modified by the mass transfer 
implicit to the proposed RLOF. If this is indeed the case, it will be 
extremely interesting to investigate the detailed surface field prop- 
erties, and to compare them to those of single O-type stars. 

A more speculative possibility is that the magnetic field was 
generated as a consequence of the RLOF event. While strong shear 
is expected to have been present during the mass transfer, it is not 
clear if the conditions required to driv e a dynamo would have been 
satisfied. In this context, we note that lTout et al.l d2008l) have pro- 
posed that highly magnetic white dwarfs have their origin in stellar 
mergers. The extension of this idea to explain the magnetic fields 
of non degenerate A, B and O type stars has been briefly described 
bv lFerrario et alj J2009I) . While the RLOF mass transfer and stellar 



Using the radius, mass and the maximum rotational period, we compute 
that the current rotation of the secondary corresponds to approximately 55% 
of the critical rate (i.e. w/oj C rit = 0.55). 



mergers are very different events, the idea that highly dynamic and 
transformative events in the evolution of A, B and O stars might be 
the origin of their magnetic fields is provocative. Given that such 
even ts are now recog nized to occur quite regularly in the Galaxy 
(e.g. lSana et alj2ol3) . these suggestions should certainly be inves- 
tigated in more detail. 

The uncertainties related to the age, distance, luminosities and 
chemical abundances of this system are a serious impediment to a 
more detailed understanding of its history and current evolution- 
ary state. Resolution of some of these issues may need to await the 
Gaia mission. In the mean time, new observations in the optical and 
X-ray domains must be acquired. In particular, higher-SNR X-ray 
spectroscopy will constrain the X-ray emitting geometry and help 
to identify the emission mechanism. Phase-resolved spectroscopy 
over the entire optical domain will allow abundances of a large 
number of chemical elements - in particular tracers of mixing and 
mass transfer - to be determined for both components. Finally, ad- 
ditional polarimetric monitoring will allow the determination of the 
rotational period of the secondary star and the detailed geometry of 
its magnetic field. 

Plaskett's star represents a heretofore unique stellar system - 
a close, massive binary with a rapidly rotating, magnetized compo- 
nent - that will be a rich system for further study. It will continue to 
be a priority target for MiMeS observations in the fall and winter of 
2012, as well as a key focus of the new Binarity and Magnetic Inter- 
actions in various Classes of Stars (BinaMiCS) project and Large 
Program at the CFHT 
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